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Abstract. An effective 517(3) x SU (3) chiral lagrangian, which includes scalar resonances, is used 
to describe the process D + — > K~n + n + at low-energies. Our main result is a set of five 5-wave 
amplitudes, suited to be used in analyses of production data. 
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INTRODUCTION 

The pioneering E791 experiment[l] has shown that heavy-meson decays are reliable 
sources of information about scalar resonances. Data about these decays are encoded 
into Dalitz plots and have to be interpreted with the help of theoretical ansatze which, 
normally, rely on Breit-Wigner expressions. In the framework of field theory, one knows 
that expressions of this kind arise naturally when two-body interaction kernels are 
unitarized. However, usual Breit-Wigner functions are problematic, since they are based 
on kernels which do not describe well the meson-meson amplitude close to threshold. 
We propose to cure this problem by means of generalized Breit-Wigner expressions, 
based on interaction kernels which comply with low-energy chiral theorems. These 
alternative trial functions can be used as tools in analyses of the decay D + — > K~n + ii + . 

Our theoretical description is based on an effective SU (3) x SU (3) chiral lagrangian, 
which includes scalar resonances and also allows a consistent treatment of the primary 
weak vertex. Final expressions represent a compromise between reliability and simplic- 
ity, so that they could be employed directly in data analyses. 

RESULTS 



The decay D + (P) — > tz + (q)iz + (q')K~ (k) involves both the weak quark transition c — ► 
sW + and strong processes associated with final state interactions. We only consider 
hadronic degrees of freedom and need the following weak vertices: (D — > nKW) a , 
(D — ► kW) a , (W — > 7t) a , (D — > KW) V , and (W — > kk) v , where k is the kappa- 
resonance and the labels v and a refer to vector and axial currents involved in W- 
couplings. Strong interactions determine most structures observed in Dalitz plots and, 
in particular, the widths of resonances. The final state one is considering contains three 
mesons and a full treatment of their interactions is impossible. Therefore we work in the 
approximation in which one of the final mesons acts as a spectator. As the emerging 



pions have isospin 2, their interactions can be neglected. Strong processes are then 
restricted to the TtK subsystem and described using the leading chiral SU (3) x SU (3) 
lagrangian given by Gasser and Leutwyler[2], complemented with resonance couplings 
from the work of Ecker, Gasser, Pich and De Rafael [3]. 
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FIGURE 1. (Color online) Diagrams contributing to the decay D + — > K n + n + ; (a) corresponds to a 
direct process, (b-e) involve the nK scattering amplitude T, and (f) depends on the production amplitude 

n. 

Our model is defined by fig. 1 . Diagram (a) represents a non-resonant background, 
since the outgoing mesons reach the detectors without interacting. The other diagrams 
involve the strong amplitudes T and EL The former describes elastic scattering and is 
unitary. The latter implements the propagation and decay of the K produced directly at 
the weak vertex[4] and has the dynamical structure shown in fig. 2. 



FIGURE 2. (Color online) Resonance propagation (full line) and decay into n + K (dashed lines); T is 
the unitary 7=1/2 scattering amplitude. 



Individual contributions from diagrams of fig. 1 to the D + decay width read 
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where G F and 0^ are the Fermi coupling constant and the Cabibbo angle, F is the 
pseudoscalar decay constant, c c \ and c m describe two possible kkK couplings, and the 
[Si] implement SU (4) symmetry breaking at the weak vertices. One has two possible kK 
subsystems and their invariant masses are jj.^ K = (q + k) 2 and fi^ , = (q' + k) 2 . In the 

construction of Dalitz plots, one first evaluates srf = [s$ a H sff] and then symmetrize 

with respect to <-> [l^ , . Allowed values for these invariant masses lie in the interval 

0.40GeV 2 < Hkx^kx' ^ 2.99GeV 2 . 

The functions Cl and T are important substructures in our results and describe respec- 
tively two-meson propagation and kK elastic scattering. The labels 1 /2 and 3 /2 refer to 
TtK isospin channels and we recall that only the former can couple with a s-channel K. 
The main features of these functions are summarized below. 
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FIGURE 3. (Color online) Bethe-Salpeter equation for the elastic kK amplitude: full equation (top) and 
perturbative solution (bottom). 



The elastic nK scattering amplitudes 7} are derived by means of the Bethe-Salpeter 
equation, represented in figs. 3, together with its perturbative solution. This geometrical 
series involves the two-body irreducible kernels and the regularized two-meson 
propagators Clj. As pointed out by Oiler and Oset[5], at low-energies, the Bethe-Salpeter 
equation can be solved exactly, and one has 

~ * (6) 



Propagators are complex above threshold and given by D.j = Rj + il. The scattering 
amplitudes are unitary and can be rewritten as 

Ti = — sind/e' 0/ tan 5/ = - - . (7) 
p l+RjJtfj 



where 8/ are real phase shifts and p = y 1 — 2 (Mj.+M 2 ) / s + (M^-M 2 ) 2 /s 2 . Our uni- 
tarization of 7} generalizes the ^-matrix approximation, which amounts to neglecting 
the real parts of Clj. 
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FIGURE 4. (Color online) Tree-level ttA" amplitude; dashed and full lines represent respectively pseu- 
doscalar mesons and scalar resonances. 

The dynamical content of the scattering amplitudes is incorporated into the kernels 
J^j, which are obtained by projecting out the 5-wave content from the diagrams of fig.4. 
Contributions from t and u channels are very small and can be neglected. Our chiral 
lagrangians then yield 

= 4^2 [4j-3p/2-4(Mj+Mi)] 
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TWO-MESON PROPAGATOR 

The propagator for a system with total momentum X = q+k, is 

_. f d 4 £ 1 
a = l Jj2^ [(l+X/2) 2 -M 2 ] l(£-X/2) 2 -M 2 ] ^ 



and can be written as Q. = — (1/167T 2 ) [L(s)+A x> ], where the function L(s) is given 
below and A,*, is an infinite constant that has to be removed by renormalization. Above 



threshold, one has 
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The renormalized propagators are given by &i = — ( 1 / 1 6 ft 2 ) [L(s) +cj] and the constants 
c/ depend on the isospin channel. They are chosen by tuning the predicted phase shifts, 
eq.(7), to experimental results at a given point s = sj. For the channel / = 1/2, the 
constant c 1/2 is chosen at the point s l/2 = m 2 K , where the experimental phase is n/2. One 
writes 

= - i^s sfcM - L ( m M ( 12 ) 

and, by construction, ^ 1/2 (m^.) = 0. In the 1 = 3/2 channel, fit to data requires c 3/2 ~ 

14.5. One notes that the functions Clj introduce other phases into the problem, given by 
tan (Oi(s) = I/Rj. 
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FIGURE 5. (Color online) Full (black) and partial contributions to the real (continuous line) and 
imaginary (dashed line) components of the amplitudes j^; the vertical scale has to be multiplied by the 
weak factor Gp cos 2 9c- 

The purpose of this work is to map the relevant degrees of freedom of the amplitude 
D + — ► K~iz + tz + at low energies. This is achieved in eqs.(l-5), which represent indi- 
vidual contributions from the diagrams shown in fig. 1 . Masses and coupling constants 



in these expressions are kept free, so that their values can be extracted from experi- 
ment. However, in order to discuss qualitative features of our results, we need to fix 
somehow these free parameters. In this case, we choose: F = 0.093 GeV, m K = 1.2 GeV, 
(c d ;c m ) = (3.2; 4.2) x 10~ 2 GeV[3] and 8 a = 8 b = 8 C = 1. We show, in fig.5, the real and 
imaginary components of the amplitudes s^i and note that the magnitudes of all contri- 
butions are comparable. Nevertheless, it is also possible to see that diagrams (a) and (f), 
which describe the non-resonant background and the direct production of the resonance 
at the weak vertex, are especially important. The latter is very sensitive to the resonance 
coupling constant cj and we take this conclusion with reserve. 

In our calculation, the non-resonating background is given by a single and very simple 
diagram. It does not contain loops and hence srf a is necessarily a real function. As the 
inclusion of loops is associated with higher orders in the chiral expansion, we conclude 
that, at low-energies, there is no phase associated with the background. It is also worth 
noting that srf a is linear in the invariant masses and its distribution over a Dalitz plot will 
not be uniform, as sometimes assumed. 

Other diagrams are complex and the presence of imaginary components can be un- 
ambiguously traced back to loops in the propagators Q/. By means of eq.(6), they blend 
with the real kernels J£/ and give rise to unitary %K scattering amplitudes. Direct in- 
spection of eqs.(2-4) shows that their phases depend on combinations of <5/ and Ofy. In 
digram (f), representing the direct production of the K-resonance at the weak vertex, the 
phase is just <5 ]/2 , as discussed in in ref.[4]. 

We assume our results for the ^ to be reliable up to values of \x^ n just above the 
point s = m 2 , corresponding to 5 l/2 = n/2. In the case of fig.5, this means invariant 

masses between 0.4 GeV 2 and 1.6 GeV 2 . Constraints imposed by chiral symmetry are 
very important at the lower end of this interval. Full details of these calculations will be 
presented elsewhere. 
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